Encouraging moderation: Clues from a simple model of ideological conflict 
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Some of the most pivotal moments in intellectual history occur when a new ideology sweeps 
through a society, supplanting an established system of beliefs in a rapid revolution of thought. Yet 
in many cases the new ideology is as extreme as the old. Why is it then that moderate positions so 
rarely prevail? Here, in the context of a simple model of opinion spreading, we test seven plausible 
strategies for deradicalizing a society and find that only one of them significantly expands the 
moderate subpopulation without risking its extinction in the process. 



The social history of ideas involves the frequent replay 
of a single story: there is a widely accepted and deeply 
ingrained dogma in the community This dogma helps to 
justify the community's institutions and shape its com- 
mon practices. Then, in the midst of this stable milieu, 
a new doctrine emerges. Backed by a small group of 
unwavering advocates, it challenges the status quo and 
steadily wins converts, eventually replacing the previous 
system to become the dominant ideology of the group. 

In some cases, there is an enduring consensus that 
the new doctrine marks a tangible improvement on the 
old. This is the case for the American civil rights move- 
ment [1], women's suffrage [2, 3], and paradigm shifts 
in science [4-6]. However in many other situations, the 
newer doctrine is not clearly better. After some time as 
the dominant approach, it too is overtaken by a younger 
alternative, which in turn is itself replaced, and so on. 
This second situation is often seen in rapidly spreading 
political campaigns [7], the booms and busts of credit 
lending and consumer confidence [8, 9], cultural fash- 
ions and short-lived reforms (e.g., Prohibition in the 
United States) [10, 11], methodological or topical fads 
in academia, and various political revolutions [12]. 

A natural question is this: why do communities, and 
sometimes entire societies, get caught in these swings 
from one ideological extreme to the other when neither 
delivers a sustainable solution? Why does not a majority 
of the population settle on an intermediate position that 
blends the best of the old and new? 

There are several ways in which this question might be 
answered, but here we give one that is purely mathemati- 
cal: the environment of successive ideological revolutions 
is not conducive to moderate mindedness simply from a 
dynamical perspective. In particular, almost all of the 
intuitive ways of encouraging moderation either fail to 
expand the moderate subpopulation or make it vulnera- 
ble to collapse in the process of encouraging its growth. 

In this Letter, we provide evidence for this claim by 
studying a minimal model of ideological revolution. Crit- 
ically, this model only addresses large-scale ideological 
conversions and does not treat the many other common 
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FIG. 1: Model structure (see text for definitions of A, B, 
AB, and A c ). Labels on the arrows indicate the allowed af- 
filiation^) of a speaker that converts a listener from one sub- 
population to another in the direction of the arrow. 



processes found in real communities, such as apparent 
conversions within the old paradigm and situations where 
there is no conversion at all but rather a splitting of opin- 
ions, or fragmentation. 

The model (Fig. 1) starts from an assumption of a 
community consisting of four nonoverlapping subpopula- 
tions: those that currently hold an extreme opinion A; 
those that currently hold the conceptually opposing opin- 
ion B (in the simplest case, just the negation of A); those 
that currently hold neither A nor B (the moderates) ; and 
those that hold A indefinitely and are immune to the in- 
fluence of others (we call these committed believers or 
A zealots). We partially overload notation, using A, B, 
AB, and A Cl respectively, to denote both the individu- 
als in these four subpopulations and the subpopulations 
themselves. This model builds on earlier pioneering work 
in sociophysics [13-20] and is directly inspired by (but 
different from) a model examined in a recent study of 
opinion dynamics [21, 22]. 

The dynamics of the basic model are determinis- 
tic, continuous and mean-field, derived as the large- 
population limit of the following random process: time is 
discrete, and at each time step we select two individuals 
uniformly at random and randomly choose one of them to 
be the speaker and the other the listener. If the speaker 
is an A or B and the listener is a B or A, respectively, 
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TABLE I: Interactions that change the membership of sub- 
populations A, B, and AB in the basic model. The A c sub- 
population is constant. 



Speaker 



Listener 
pre-interaction 



Listener 
post-interaction 



A, A c 



B 



AB 



AB 



A 



B 



A 



AB 



AB 



B 



then the listener is dissuaded from his or her extremist 
position and becomes an AB. However, if the listener is 
an AB, then the listener becomes an A if the speaker is 
an A and a B if the speaker is a B. In all other cases, 
there is no change in the state of the speaker and listener 
(Table I). Note that in this highly simplified model, mod- 
erate speakers do not produce a change of mind in cither 
their listeners or themselves; only extremists successfully 
rally others to their cause. 

Let Ha, n>B, and uab denote the expected fractions 
of the total population of N individuals corresponding 
to the uncommitted A, B. and AB subpopulations, re- 
spectively, and let p denote the constant fraction of the 
population in the committed A c subpopulation. We will 
study how varying p, the proportion of zealots, affects 
the eventual state of the system. Using this new nota- 
tion, we can consider the expected change to the subpop- 
ulation fractions in the limits of a large population and 
a vanishing time step (which we take to grow like N and 
shrink like 1/iV, respectively). This reduces our discrete 
dynamics to the following rate equations: 



flA = (j> + n A )n A B - n A n B , 
riB = n B n AB - (p + n A )n B , 



(1) 



where uab = 1 — P~ tia— n B and the overdot denotes dif- 
ferentiation by time. Since we present no formal evidence 
that the dynamics of (1) do actually occur in practice, our 
work could alternatively be viewed as posing this model 
and its subsequent generalizations as interesting in their 
own right. 

Now suppose we run the system (1) to equilibrium 
starting from a population initially composed of only A c 
and B individuals. We will use this initial condition for 
all the systems considered in this Letter; the idea is that 
A represents the new doctrine and B the reigning view. 
If we then track the final fractions of ua, n B , and n A B as 
functions of p, we find (as in Ref. [21]) that the equilib- 
rium state changes dramatically as we increase p through 
a critical value p c (Fig. 2) . For p < p c the system remains 
similar to how it started — most of the individuals main- 
tain B. However as p is increased through p c , the system 
undergoes a discontinuous transition, and for p > p c the 
entire population quickly reaches a consensus on A. A bi- 
furcation analysis shows that p c — 1 — V3/2 ss 0.134 [23]. 




FIG. 2: The equilibrium values of ua (red plus signs), nj 
(blue dots), and uab (magenta open circles) for the basic 
model as functions of p, assuming an initial population with 
(nA,ne) = (0, 1 — p). The vertical dashed line marks the 
critical value p c = 1 — v3/2 ~ 0.134. At values of p greater 
than p c , riB, and uab are zero and ua = 1 — p- 



To test the robustness of these mean-field predictions, 
we simulate the model on a diverse set of real social net- 
works. Figure 3 shows that in each case, the n B vs p 
curves resemble the mean-field result depicted in Fig. 2. 
The primary differences are a lower p c value for the real 
networks and a small, stable fraction of peripherally lo- 
cated B believers for p > p c . 

With (1) as our starting point, we now ask how we 
might alter the model to encourage moderation. Specifi- 
cally, we would like to (i) increase the equilibrium size of 
the moderate subpopulation, and (ii) decrease the chance 
that this equilibrium size could drop substantially if the 
parameter values (just p for the basic model) were to vary 
a little. In search of a strategy that does both, we explore 
seven different generalizations of the basic model. Three 
generalizations are discussed here in the main text; the 
rest are treated in the Supplemental Material [23]. Fig- 
ures 4(a)-(c) summarize mean-field results for these three 
generalizations, and for comparison, the corresponding 
simulation results on a real social network arc shown in 
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FIG. 3: The equilibrium fraction of B believers remaining af- 
ter the basic model is run on the giant connected component 
of (a) the U.S. network of corporate board memberships in 
1994 [24], (b) four coauthorship networks in the physics di- 
vision of arXiv.org [25], and (c) the friendship networks of 
the location-based social networking websites Gowalla and 
Brightkite [26]. Note the abrupt transitions in ns; compare 
with the corresponding curve in Fig. 2. 
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FIG. 4: (a)-(c): Mean-field results obtained analytically for generalizations of the basic model with (a) stubborn moderates, (b) 
evangelical moderates, and (c) nonsocial deradicalization (see the text for details). The final equilibrium values of n A (red plus 
signs), tie (blue dots), and uab (magenta open circles) for the initial condition (riA,ns) = (0, 1 — p) are plotted as a function 
of the new parameter (s, r, or u) in the corresponding generalized model. Of the three strategies shown — and in fact for all 
seven considered in the Supplemental Material [23] — only nonsocial deradicalization allows for the growth of the moderate 
fraction up to 1 — p without risking its extinction, (d)-(f): Representative simulation results for the discrete-time versions of 
the models with (d) stubborn moderates, (e) evangelical moderates, and (f) nonsocial deradicalization when these models are 
run on the arXiv coauthorship network for high energy physics theory (hep-th) [25]. The plots show the equilibrium fractions 
of ua (red plus signs), ub (blue dots), and uab (magenta open circles) obtained. Each simulation is started from the state in 
which a random but highly interconnected fraction p of the population is committed to a belief in A and the rest believe B. 
The simulation is then run for 10 8 time steps after which the values of ua, riB, and uab axe tabulated. The constant fractions 
of zealots in the six panels of this figure are (a) 0.1, (b) 0.1, (c) 0.05, (d) 0.035, (e) 0.02, and (f) 0.02. 



the panels beneath them [Figs. 4(d)-(f)]. Importantly, 
these do not constitute full empirical validations of the 
model and its generalizations (which would require dy- 
namical data that are hard to obtain). Rather, we in- 
clude these simulations only as an indication of where 
the results of such tests might lie. Furthermore, we only 
consider the equilibrium values reached from the pre- 
rcvolution initial condition (n^,^) = (0, 1— p)\ as dy- 
namical systems, the basic model and its generalizations 
are capable of a wider range of behaviors [23] . 

As a first attempt at achieving (i) and (ii) above, sup- 
pose we could somehow make the moderates less likely 
to convert to either of the two radical positions. We can 
represent this by generalizing the basic model to 

n A = {I - s)(p + n A )n AB -n A n Bl 
n B = (1 - s)n B n AB - (p+ n A )n B , 

where the stubbornness parameter s indicates how likely 
a moderate is to remain moderate after listening to an 
extremist. When s is zero, we recover the basic model. 

Intuitively, one might expect that increasing s should 
increase the size of the moderate subpopulation. Indeed, 
when s is small enough, the moderate subpopulation does 
grow slightly with increasing s [Figs. 4(a) and 4(d)]. But 



remarkably, if s increases past a certain threshold, the 
moderates are driven to extinction; the size of their sub- 
population drops to zero. 

We can examine this surprising behavior in another 
way by calculating how s affects p c (the critical fraction 
of zealots needed for the revolution to succeed). Intuition 
would suggest that p c should increase with s; the more 
stubborn the moderates are, the more zealots are needed 
to persuade them and everyone else. In fact the opposite 
is true: p c decreases with s, dropping monotonically from 
1 — V3/2 at s = to zero at s = 1 [23]. Thus, increasing 
the stubbornness of the moderates makes the population 
more vulnerable to takeover by the zealots. 

To make sense of why p c should decrease with s, it 
helps to realize that increasing s not only reduces the 
flow of AB individuals to opinion A but also to opin- 
ion B, thereby depleting both the uncommitted A and 
B subpopulations. With competition from B extrem- 
ists over the AB subpopulation weakened as a result, it 
takes fewer A zealots (and hence a lower p c ) to convert 
the moderates to the A camp. 

This explanation suggests that evangelism is an impor- 
tant force in the dynamics. So as a second strategy, we 
might try having the moderates actively promote mod- 
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eration via the following generalization: 

flA = (p + riAjTLAB - TlATlB - mATlAB, 

hs = n B n A B - (p + n A )n B - rn B n A B, 



(3) 



where the new parameter r is a nonnegative real number 
that reflects the intensity of the moderates' evangelism. 

Again it may seem intuitively clear that the size of 
the moderate fraction should increase if the moderates 
start actively deradicalizing the population. For r up to 
unity, however, the outcome is similar to that of mak- 
ing the moderates more stubborn. Figures 4(b) and 4(e) 
show that at a certain value of r, the size of the moder- 
ate subpopulation snaps discontinuously to zero. If the 
moderates' campaign of persuasion is sufficiently success- 
ful from the start — that is, if r starts and stays large 
enough — then the moderates do in fact maintain a large, 
robust equilibrium population. However, if they fail 
to sustain this level of persuasiveness indefinitely, their 
evangelistic efforts can instigate their own extinction. 

Finally, let us consider a third strategy: suppose that 
the fanatics are deradicalized by a promoderation media 
campaign or other environmental stimulus rather than 
through social interaction with moderates. We could 
then expect the dynamics to take the following form: 



h.A = (p + n A )n A B - n A n B - un A , 
h B = n B n AB ~{p + n A )n B - un B , 



(4) 



where it is a nonnegative parameter representing the rate 
at which the radicals abandon their radical position in 
response to the nonsocial stimulus. 

In contrast to the first two strategies (as well as four 
others treated in the Supplemental Material [23]), in- 
creasing the new parameter (u) in this system generally 
increases the equilibrium uab toward a limit of 1 — p, 
with the one exception of a discontinuous drop partway 
through the ascent in Fig. 4(c). However, the drop is 
not to zero as it was for the other strategies, and it van- 
ishes in the limit of small p. Furthermore, following the 
drop, regrowth of uab is rapid. Hence, this mechanism 
of promoting moderation, which we might call nonsocial 
deradicalization, provides the first acceptable means that 
we have found for expanding the moderate population 
in the midst of an ideological revolution. This holds for 
the three strategics in this Letter, and also for the four 
others in the Supplemental Material [23]. 

By itself, this final assessment should be regarded with 
caution. We suggest that a greater emphasis be placed on 
our general approach as a framework for testing possible 
strategies as part of a continuing research program, which 
through further study might well uncover other means 
of fostering moderation more sophisticated than those 
considered here. 
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Supplemental Material 



Here we derive the mean-field results described in the 
main text. Our approach consists of analyzing seven dif- 
ferent one-parameter generalizations of the basic model. 
Our central conclusion is that only one of these, nonsocial 
deradicalization, offers a robust means of expanding the 
moderate subpopulation. We finish with an example of a 
two-parameter generalization that, like nonsocial derad- 
icalization, can deradicalize all uncommitted extremists 
without risking the extinction of the moderates. 

In order to help interpret the results below, we recom- 
mend that the reader plot the nullclines of each system 
and study how they vary with the parameter values using 
the MATLAB program pplane. This program is avail- 
able online at http://math.rice.edu/~dfield/dfpp.html 
and is free for educational use. 



THE BASIC MODEL 

The dynamics of the basic model are 

h A = (p + n A )n AB - n A n Bl 

h B = n B n A B — (p + n A )n Bl 

where the overdot represents differentiation by time and 
n AB = I — p — n A — n B (see the main text for definitions 
of the remaining terms). The nullclines of (1) are 

= (p + n A )n A B - n A n B , (2a) 
= n B n AB - (p + n A )n B . (2b) 

The first implies that 

{p + n A )(l - p- n A ) 



(1) 



n B 



(3) 



p + 2n A 

and the second (by factoring) implies either n B = or 
n B = l-2p- 2n A . (4) 

If we substitute n B = into (2a), we obtain two solu- 
tions: n A = —p and n A = 1 — p. So (f ) has fixed points 
at (n A ,n B ) = (—p,0) and (n A ,n B ) = 0). 

To obtain the remaining fixed points, we set (3) and 
(4) equal to each other and solve the resulting quadratic 
in n A to obtain 



Ul-4p± VV-8p+f 



(5) 



Together with (4), this gives a third and fourth fixed 
points so long as the discriminant of (5) is nonnegativc. 

Solving for the roots of the discriminant of (5) gives 
p c = 1 ± V3/2. We can disregard the value with the 
positive sign because p represents a fraction of the total 
population and thus must be less than one. The value 
with the negative sign represents the p at which the third 
and fourth fixed points coalesce in a saddle-node bifur- 
cation in (n Al n B ) phase space. 



GENERALIZATION I: 
STUBBORN MODERATES 

One way to encourage expansion of the moderate sub- 
population would be to try to make existing moderates 
less likely to abandon their moderate views and convert 
to an extreme position. This idea can be captured with 
the following generalization of (I): 



h A = (I - s)(p + n A )n AB - n A n B , 
h B = (1 - s)n B n AB - (p + n A )n B . 



(6) 



Here s is a parameter on the unit interval. Increases in 
s correspond to the moderates becoming less prone to 
switch to a radical position. 

Following the approach of the previous section, we find 
that this system has fixed points at (n Al n B ) ~ (— p, 0) 
and (n Al n B ) = (1 — p, 0), as well as where 



(l-s)(l-3p)-p±y/A\%p) 

UA = R A 7 ) 

6 — 4s 



and 



(I -. s)(l -p)-p- (2 -s)n A 
n B = — s ' (8) 

with A(s,p) = (2 - s)V - 2(1 - s)(4- 3s)p + (f - s) 2 . 
This second pair of fixed points exist everywhere that 
A(s,p) > 0, and the curve A(s,p) — represents the 
parameter pairs (s,p) at which the third and fourth fixed 
points coalesce in a saddle-node bifurcation. 

Solving A(s,p) = in terms of p (i.e. finding the roots 
of the quadratic polynomial in p) gives 



Pc 



(f - s)(4 - 3s) ± 2^(f-s) 3 (3-2s) 



(9) 



(The c subscript is added to indicate a connection to 
literature on critical phenomena and phase transitions.) 
The root in (9) with the plus sign represents a saddle- 
node bifurcation that occurs in the third quadrant of the 
(n A ,n B ) plane where n A and n B are both negative, so 
the root with the minus sign is the only one of interest to 
us. This curve descends monotonically fromp = 1— v3/2 
at s = to p = at s = f . For the (s,p) pairs below it, 
A(s,p) > 0, so the third and fourth fixed points exist for 
these (s,p). The third and fourth fixed points generally 
have small values for n A and large values for n B , and 
the one having the n A with the negative sign in (7) is 
stable. (We can show this by an analysis of the Jacobian 
of (6).) By contrast for the (s,p) pairs above the curve, 
A(s,p) < 0, so the only fixed point of relevance in this 
region of parameter space is (n A ,n B ) = (I — p, 0). Again 
by stability analysis on the Jacobian of (6), this point is 
stable for all (s,p) pairs on the unit square. 
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GENERALIZATION II: 
EVANGELICAL MODERATES 



GENERALIZATION III: 
NONSOCIAL DERADICALIZATION 



Another approach to increasing the moderate fraction 
would be to make the moderates evangelical for their 
moderate perspective. To reflect this idea, we can gener- 
alize the basic model in the following way: 

h A = (p + n A )n AB - n A n B - rn A n ABl 



fiB = n B n AB - (p + n A )n B - rn B n AB , 



(10) 



where r is a nonnegative parameter representing the ef- 
fectiveness of the moderates' evangelism. 

Analyzing this model as we did the basic model, we 
again find that there are two fixed points on the n A axis: 
(n A ,n B ) = (1-p, 0) and (n A ,n B ) = (-p/(l-r),0), and 
two fixed points off the n A axis with locations given by 



(1 - r)(l - 3p) - (1 + r)p ±VA^P) , in 



6-4r 



and 



n B = 



(1 - r)(l - p) - p - (2 - r)n A 
1-r 



(12) 



where A(r,p) = 4(1 - r) 2 p 2 - 4(1 - r)(2 - r)p + (1 - r) 2 . 

The curve A(r,p) = gives the (r,p) pairs at which 
a saddle-node bifurcation of the third and fourth fixed 
points occurs. However since p must be on the unit in- 
terval and r is nonnegative, it can be easily shown that 
the only relevant segment of this curve is 



Pc 



2-r- v/3 - 2r 



(13) 



2 - 2r 

where r is restricted to the unit interval. This segment 
starts at p = 1 — >/3/2 when r = and declines mono- 
tonically to zero as r reaches 1. For allowed (r,p) below 
the curve, A(r,p) > 0. So the two fixed points given 
by (11) and (12) exist in this region. For allowed (r,p) 
above the curve (13), A(r,p) < and so (11) and (12) 
do not exist in this region and the only fixed points are 
(n A ,n B ) = (1 -p,0) and (n A ,n B ) = (-p/(l - r),0). 

There is also a second bifurcation, this time transcrit- 
ical, that occurs when the fixed point on the n A axis 
located at n A = —p/(l — r) passes leftward through 
the fixed point at n A = 1 — p as r is increased through 
r = 1/(1 — p). For comparison with (13), we can write 
this bifurcation curve as a critical p: p c = 1 — 1/r. This 
curve starts from p = at r — 1 and asymptotically ap- 
proaches p = 1 as r is increased toward positive infinity. 

From analysis of the Jacobian of (10), we find that 
the fixed point with the n A in (11) having the negative 
sign is stable, while the other is unstable. We also find 
that the fixed point at (n Al n B ) = (1 — p, 0) is stable 
when < r < 1/(1 — p). However as r is increased 
through the transcritical bifurcation, the fixed point at 
n A = —p/(l — r) inherits the stability of the fixed point 
at 1 —p. Thus for r > 1/(1 —p), the equilibrium value of 
n A decreases to zero like 1/r for large r. 



A third way to try to increase the size and stability of 
the moderate subpopulation would be to provoke radi- 
cals on both sides to become moderate by some nonso- 
cial means. This strategy is captured by the following 
generalization of the basic model: 

fiA = (p + n A )n AB - n A n B - un A , 
hs = n B n AB - (p + n A )n B - un Bl 

where u is a nonnegative parameter representing the rate 
of dcradicalization. 

We can again treat this system in the same way that 
we did the basic model. The first pair of fixed points 
have n B = and 



(14) 



U = \ (l - 2p - u ± \/4up+(l-M) 2 ) . 



(15) 



Note that the discriminant of (15) is always positive, so 
these first two fixed points exist for all allowed p and 
u. For fixed p, the positive root of (15) shrinks to zero 
like 1/u as u is increased toward positive infinity. This 
is easily seen from the asymptotic expansion of n A in 
powers of 1/u, which begins 



.1 (1 - 2p) (1 - 2p) 2 

n A =p(l-p) 



(16) 



For the second pair of fixed points, the first coordinate 
is given by 



n A 



i ( 1 - Ap - u ± y/A(u,p) 



(17) 



and the second coordinate by 

n B = l-2p-u-2n A , (18) 

where A(u,p) = 4p 2 - 4(2 + u)p + (1 - u) 2 . 

Solving A(u, p) = for p, we obtain the curve on which 
the saddle-node bifurcations occur: 



Pc 



1 



u v3 + 6u 
2 2 ■ 



(19) 



Only the root in (19) with the minus sign is of interest 
to us, because only this root has values less than 1 for 
positive it. Additionally, we are only interested in the 
segment of this curve with u on the unit interval, since 
for larger u the curve represents a saddle-node bifurcation 
in the third quadrant of the (n A , n B ) plane. 

Below this curve segment, we find that A(u,p) > 0, so 
the third and fourth fixed points exist there. By analyz- 
ing the Jacobian of (14), we can show that the fixed point 
having the n A in (14) with the negative sign is stable and 
draws in trajectories to a small n A and large n B . How- 
ever above the curve segment, we find that A(u,p) < 0, 
so the only fixed points of (14) are the two with n B = 
and n A values given by (15). The one of these with the 
positive n A value is stable and draws in all trajectories 
of interest for all allowed (u,p) above the curve. 
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GENERALIZATIONS IV, V, AND VI: 
ASYMMETRIC VERSIONS OF 
GENERALIZATIONS I, II, AND III 

Generalizations I and II above lead to an annihilation 
of the moderate subpopulation for some values of s and 
r. This occurs because at these values of s and r, much 
of the B subpopulation is absorbed into the moderate 
subpopulation early in the time series, which allows the 
A believers to subsequently convert the entire popula- 
tion to their belief A. Motivated by this observation, we 
now consider three variations on the above generaliza- 
tions that all attempt to shift the dynamics in favor of 
the B believers. We first show that each of these models 
has a fixed point that converges to (jia, tib) = (0, 1 — 2p) 
in a large-parameter limit. We then discuss the implica- 
tions of this for the fraction of moderates in each model. 

To prepare the first variation (Generalization IV), we 
omit the factor of (1 — s) in the second equation of (6) 
to obtain 



fiA = (1 - s)(p + n A )n A B - n A n B , 
h B = n B n A B — (p + n A )n B - 



(20) 



Like the original system, this new system has fixed points 
at (nA,ns) = (— p, 0) and (n^ng) = (1 —p,0). It also 
has fixed points where 



1 - 4p + 2sp ± ^4p 2 - 4(2 -s)p + l 
nA = ^5 21) 

6 — 2s 

and tib = 1 — 2p — In a- In the limit as s tends to 1, this 
second pair of points converges to (ua, Ub) = (1/2 — p,0) 
and (n A , n B ) = (0, 1 - 2p). 

A second variation (Generalization V) is obtained by 
discarding the last term in the second equation of (10). 
This yields 



n A = (p + n A )n A B - n A n B - rn A n A B, 
Ub = ubUab - (p + n A )n B - 



(22) 



The fixed points of this system are located at (ua, tib) 
{—p/ (1 — r), 0) and (uajTIb) = (1 — P, 0), and where 



1 - Ap + rp ± y/(2 - rfp 2 - 2(4 - r)p + 1 

n A = 7. 7, (23) 

6 — 2s 

and ub = I — 2p — 2ua- In the large-r limit, the fixed 
point having the ua with the minus sign in (23) becomes 
(n A ,n B ) = (0,l-2p). 

Finally, the third variation (Generalization VI) in- 
volves omitting the final term in the second equation of 
(14) to obtain 



n A = {p + n A )n A B ~ n A n B - un Al 



h B = n B n A B ~ (P + n A )n B - 
This system has two fixed points along the nA axis at 



(24) 



It also has an additional two where 
n A = I (l-4p+u± VV - 8(1 + u)p + (1 + u) 2 ) , (26) 

and nB = l — 2p—2nA- Analogous to the two cases above, 
the fixed point with the negative sign in (26) converges 
to (nA,n-B) = (0, 1 — 2p) in the large-w limit. 

The point (n^ns) = (0, 1 — 2p) therefore appears as 
the limit of a fixed point in all three of these new models. 
Furthermore at this fixed point, the moderate subpopu- 
lation is nAB = P, so neither p nor nAB can exceed 1/2 
at this fixed point. A more extensive analysis of the ex- 
istence and stability of all fixed points shows that when 
p < 1/2 and each parameter (s, r, and u) is large enough, 
this fixed point exists in the first quadrant and is stable, 
receiving the trajectories of a large basin of attraction 
around it including the initial point (ua, tib) = (0, 1 —p). 



GENERALIZATION VII: 
OPPOSING ZEALOTS 

A seventh strategy for expanding the fraction of mod- 
erates in the population is to assemble a faction commit- 
ted to B that could counterbalance the existing faction 
committed to A. If we represent the size of this second 
fraction by q, then the basic model becomes 



h A = (p + n A )n A B - (q + n B )n A , 
n B = (q + n B )n A B - (p + n A )n B , 



(27) 



n A 



= 5(1 - 2p- u± V / Aup+ (1 - u) 2 



(25) 



where now nAB = 1 — P — 1 — nA — tib ■ 

We can find exact expressions for the fixed points of 
this system, but they are more complicated than those 
of the first three generalizations above. This is primarily 
due to the fact that we cannot factor out an ns from the 
left-hand side of the second rate equation in (27), so we 
must use the quartic equation (rather than the quadratic 
equation) to obtain the fixed points. Due to the size of 
the solutions, we don't include them here. However they 
can be easily found with a computational platform for 
symbolic algebra such as Mathematica. 

When we arbitrarily vary p and q, the fixed points of 
(27) generally undergo saddle-node bifurcations in the 
(jiA^ns) phase plane. In the special case that we set 
q = p, varying the combined parameter (which we will 
just call p from now on) produces a two-dimensional su- 
percritical pitchfork bifurcation. The handle and middle 
tine of this pitchfork lie along the line n B = nA, with 
the origin located along the handle. At some critical p, 
the two outside tines of the pitchfork emerge and grow 
outward symmetrically about ns = tia- 

To find the critical p where this bifurcation occurs, 
let's consider the nullclines of (14) under the substitution 
q = p: 

= (p + n A )n A B - (p + n B )n A , (28a) 
= {p + n B )n A B - (p + n A )n B - (28b) 
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If we subtract each side of (28b) from the corresponding 
side of (28a), we obtain = (ua — tib)(1 — 3p — tia — ub), 
so either ns — Ua or ua + tib = 1 — 3p. The former, 
by substitution into (28a), gives the fixed points on the 
line n B = n A : (n A ,n B ) = {—p,—p) and (n A ,n B ) = 
((1 - 2p)/3, (1 - 2p)/3). The latter, also by substitution 
into (28a), yields = 



- (1 - 3p)n A +p 2 , or 



n A 



(29) 



Since both p and 1 — 3p must be positive (recall that 
1 — 3p is equal to ua + Ub), the discriminant vanishes 
when 1 — 3p = 2p, or when p = 1/5. So this is the p of 
the pitchfork bifurcation that we seek. 

In fact we not only can find the p of the pitchfork bifur- 
cation but also the (q, p) pairs at which the saddle- node 
bifurcations occur. We can do so by identifying the parts 
of the fixed points of (27) that are under even-powered 
roots. In fact, there is only one such part (although it 
appears many times throughout the expressions for the 
fixed points) and setting it equal to zero gives an implicit 
function for the curve in parameter space along which the 
bifurcations occur. With some simplification, this is 



32(p 3 c? + pq 3 ) + 61pV " W0{p 2 q + pq 2 ) 

+ 82pq + 4(p 2 + q 2 ) - 8(p + q) + 1 = 0. 



(30) 



Since the highest power of p and q in (30) is 3, this level 
set could instead be expressed as a set of two explicit 
functions, either p in terms of q or the other way around. 
The two separate functions correspond to two different 
saddle-node bifurcations, and the locus where they meet 
marks the p of the pitchfork bifurcation. 

More broadly, (30) represents a cusp catastrophe and 
can be shifted and rotated into a more conventional ori- 
entation (cusp positioned at the origin, pointing to the 
right) with the new parameters a = p + q — 2/5 and 
b = p — q. The variable coordinate, orthogonal to the pa- 
rameter plane, may be taken to be ua, n-B, or some rea- 
sonable combination of these, for example (ua — "■b)/v / 2 
(the projection onto the line = ua + ns)- 

A variety of interesting results can be proven for the 
system (27). For example, one is that the hab frac- 
tion at equilibrium cannot exceed (3 — V3)/6 (< 0.212). 
Showing this consists of demonstrating that in either of 
the two regions of parameter space partitioned by the 
curve (30), the fraction of uab at equilibrium is great- 
est along this curve itself, with the maximum for the 
region with small p and q occurring at the intersection of 
(30) with the p axis, that is (q,p) = (0, 1 — \/3/2), and 
the maximum for the complementary region of parame- 
ter space occurring at the locus of the cusp bifurcation 
where {q,p) = (1/5,1/5). (Here we have assumed that 
the system is started from (nA,ns) = (0, 1 — p), but the 
result holds for any initial condition.) Evaluation of the 
function for uab at these points demonstrates that the 
former has the larger value at (3 — \/3)/6. 



GENERALIZATION VIII: 
STUBBORN MODERATES 
AND OPPOSING ZEALOTS 

Finally, we analyze a model that combines Generaliza- 
tions I and VII: 



tia = (1 - s)(p + n A )n A B ~(q + n B )n Al 
h B = (1 - s)(q + n B )n A B ~{p + n A )n Bl 



(31) 



where uab is again 1 — p — q — tia — nB ■ 

Our aim in studying this model is to evaluate whether 
we can expand the moderate subpopulation with the fol- 
lowing two-part strategy: (i) assemble a faction commit- 
ted to B of the same size as that committed to A, and 
(ii) increase s (the stubbornness of the moderates to re- 
main moderate) . For simplicity, we suppose that we have 
already rallied a faction committed to B equal to the size 
of that committed to A, or, in the parlance of the model, 
q = p. 

This leaves us with the following nullclines: 

= (! - s)(p + n A )n AB - (p + n B )n A , (32a) 
= (1 - s)(p + n B )n A B - (p + n A )n B - (32b) 

We expect a pitchfork bifurcation to occur at some value 
p as a function of s. For s = 0, this p is just 1/5, as we 
found in our analysis of Generalization VII. 

If we subtract each side of (32b) from the corre- 
sponding side of (32a), we obtain = (ua — ™b)[(1 — 
s)(l — 2p — ha — tlb) — p], so either tlb = nA or 
ua + ns = 1 — 2p — p/(l — s). The former combined 
with (32a) gives the fixed points nA = ns = ~p and 
riA = nB = (1 — s)(l — 2p)/(3 — 2s), while the latter 
combined with same equation yields = n 2 A — Xn A + p 2 , 
where we have defined A 
have 



1 - 2p - p/(l - s). Thus we 



n A 



i(A±VA 2 - V). 



(33) 



Now p must be positive, and since A is the sum Ha + nB, 
it also must be positive. So the discriminant vanishes 
when A = 2p, or at 



Pc 



1-S 

5 -4s' 



(34) 



which is the p at which the pitchfork bifurcation occurs. 
Note also that at this p c we have the beautiful result that 
n A = n B = P- 

Since (34) starts at p — 1/5 when s = and descends 
monotonically to p = at s = 1, we can see immediately 
that no matter how small p is, if we increase s enough 
we will cross through the p c given by (34) and obtain an 
iiab of (1 — 2p)/(3 — 2s). Moreover, before crossing this 
Pc nAB is of size p/(l — s). With additional work, we 
can show that this implies the rise in nAB is continuous 
and therefore robust to minor variations in s and p. 
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BIFURCATION SETS AND PHASE PORTRAITS FOR MODELS I, II, III, and VII 



Model I: 
Stubborn Moderates 



Model II: 
Evangelical Moderates 



Model III: 
Nonsocial Deradicalization 



Model VII: 
Opposing Zealots 




Pc 



(0) 



(d) 
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FIG. 1: (color online) (a)-(d): Bifurcation sets for the models I, II, III, and VII with representative phase portraits shown in 
miniature. The solid brown curves indicate boundaries along which bifurcations occur. Specifically, the curves in (a) and (c), 
each branch of the curve in (d), and the left curve in (b) mark parameter pairs at which saddle-node bifurcations occur. The 
right curve in (b) signifies a transcritical bifurcation, and the cusp of the curve in (d) coincides with a supercritical pitchfork 
bifurcation. More broadly, the curve in (d) reflects a cusp catastrophe. The qualitative structure of the phase plane changes 
across these bifurcation curves. However it remains the same within each connected region and is illustrated by the inset 
phase portraits. Each phase portrait has horizontal and vertical axes of ua and ns respectively and is restricted to the unit 
square. The solid lines are nullclines and the closed and open dots represent, respectively, stable and unstable fixed points. The 
horizontal dotted lines show the trajectories followed by corresponding plots (e)-(h). From left to right, the p values of these 
dotted lines are 0.1, 0.1, 0.05, and 0.15. (e)-(h): The final equilibrium values of ua (red plus signs), tib (blue dots) and uab 
(magenta open circles) for the initial condition (n.4,ns) = (0, 1 — p) as a function of the new parameters (s, r, u, and q). Note 
that only nonsocial deradicalization allows for arbitrary growth of the moderate population without risking its extinction. All 
curves in panels (a)-(h) are computed exactly in the above sections. 



